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纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算法 
付 鹏 斌 ， 铁 惠 杰 ， 杨 惠 荣 


(北京 工业 大 学 信息 学 部 ,北京 100124) 


摘 要 : 针对 图 像 重建 过 程 中 产生 的 边缘 结构 被 破坏 和 纹理 细节 丢失 问题 ， 提 出 一 种 纹理 细节 和 边缘 结构 保持 的 图 
像 插 值 算法 。 首 先 ， 采 用 自 适 应 赋值 的 八方 向 边缘 检测 划分 图 像 区 域 ， 其 次 ， 构 造 双 变 量 有 理 吕 数 模型 ， 它 可 在 有 
理 模 型 和 多 项 式 模型 之 间 转 换 ; 最 后 提出 基于 边缘 数据 的 局 部 不 对 称 性 和 梯度 特征 调整 待 酝 值 点 空间 距离 的 方法 ， 
调整 边缘 部 分 待 插 值 点 坐标 并 代入 有 理 模型 实现 插值 ， 而 非 边缘 部 分 采用 多 项 式 模型 插值 。 实 验 表 明 ， 该 算法 的 峰 
值 信 骂 比 平均 提升 了 0.48dB-2.17d4B， 结 构 相 似 性 平均 提升 了 0.004-0.028， 获 得 了 较 高 的 客观 评价 数据 。 该 算法 将 原 
空间 距离 不 变 的 插值 修正 为 空间 距离 变化 的 插值 ， 有 效 地 保持 了 图 像 的 边缘 结构 和 纹理 细节 ， 使 得 重建 结果 具有 较 
好 的 视觉 效果 。 
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Texture details and edge structure are maintained by image interpolation algorithm 


We Fu Pengbin, Tie Huijie, Yang Huirong 
™y (Faculty of Information Technology, Beijing University of Technology, Beijing 100124, China) 


Abstract: In order to solve the problem that in the process of image reconstruction destroys the edge structure and losses the 
texture detail , this paper proposed an image interpolation algorithm for preserving the texture detail and the edge structure. 
Firstly, image divides area by using eight direction edge detection based on adaptive threshold. Secondly, construct a bivariate 
rational function model, which can convert into rational model and polynomial model. Finally, based on the local asymmetry 
edge data and gradient feature, this paper proposed a method to adjust the distance of points to be interpolated. By this method, 
adjusting the coordinates of the points to be interpolated and using rational model to realize interpolation, as for the non edge 
part, adopt polynomial model to interpolation. Experiments show that, the peak signal to noise ratio of the algorithm increases 
0.48db-2.17db on average, and the structural similarity index increases 0.004-0.028 on average. This algorithm obtains high 
objective evaluation data. In this algorithm, the original spatial distance invariant interpolation is modified to the spatial 


distance varying interpolation, which effectively keeps the edge structure and texture details of the image, and makes the 
reconstruction result have better visual effect. 
Key words: image interpolation; texture details and edge structure remain; adaptive gradient; rational function; spatial 


distance correction 


be 图 像 插值 技术 是 指 低 分 辨 率 图 像 经 过 插值 变换 转换 为 高 。 边缘 结构 的 保持 并 不 理想 。 
一 。 分辨 率 图 像 的 一 种 图 像 处 理 技术 ， 在 遥感 成 像 、 医 学 影像 、 为 实现 图 像 边缘 的 平滑 过 渡 和 纹理 细节 保 真 ， 本 文 借鉴 
公共 安全 、 计 算 机 视觉 、 数 字 娱 乐 及 多 媒体 通信 等 领域 有 着 ” 双 变 量 有 理 插值 模型 BH， 设计 了 纹理 细节 和 边缘 结构 保持 
广泛 的 应 用 。 的 图 像 插 值 算法 。 在 双 变量 有 理 插值 模型 的 基础 上 ， 基 于 边 
常见 的 插值 算法 依据 插值 方式 的 不 同 可 以 被 分 成 两 类 ， 缘 数据 的 局 部 不 对 称 性 调整 边缘 区 域 待 插值 点 空间 距离 


和 


加 灿 
过 
可 
亏 六 过 过 


离散 插值 算法 和 连续 插值 算法 。 离 散 插值 算法 是 指 根据 已 知 。” 步 确定 待 插值 点 坐标 。 随 后 通过 获取 待 插值 点 水 平和 重 
点 像素 值 之 间 的 内 在 联系 ， 建 立 某 种 变换 来 确定 未 知 点 的 像 ”向 的 梯度 角 ， 计 算 已 知 点 在 待 插值 点 水 平和 垂直 梯度 六 
素 值 。 如 常见 的 最 近邻 、 双 线性 中 等 插值 算法 , 但 是 这 类 算法 。 素 值 的 梯度 所 占 权 重 。 并 根据 获得 的 权重 实现 待 插 
的 插值 结果 会 出 现 模糊 失真 等 状况 。 连 续 插值 算法 指 的 是 应 梯度 的 空间 距离 修正 ， 获 得 边缘 区 域 待 插值 点 修 


I 
I 


T 
Tm; 
rT 
Pa 
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已 知 点 像素 值 之 间 的 联系 构造 插值 曲面 ， 计 算 未 知 点 的 像素 ” 标 。 纹 理 细 节 和 边缘 结构 保持 的 图 像 插 值 算法 ， 在 保留 原 扣 


值 。 如 常见 的 双 三 次 外、 三 次 样 条 中 等 算法 ,但 是 这 类 算法 也 ”信和 模 型 纹理 细节 保 真 度 高 的 基础 上 ， 将 原 空间 距离 不 变 的 插 
会 使 获得 的 图 像 出 现 锯齿， 模糊 等 状况 。 为 解决 已 有 插值 算 。“” 值 算法 改进 为 空间 距离 变化 的 插值 算法 ， 提 升 了 边缘 区 域 的 
法 中 存在 的 问题 ， 文 献 [4] 提 出 一 种 边缘 定向 的 插值 算法 ， 该 。 平滑 性 。 实 验 结果 表明 ， 本 文 所 提 算 法 的 重建 结果 不 仅 获得 
方法 具有 较 好 的 边缘 重建 效果 ， 但 在 纹理 细节 部 分 容易 产生 。 了 最 高 的 峰值 信 噪 比 均值 和 结构 相似 性 均值 ， 具 有 较 高 的 客 
走样 、 噪 声 问题 。 为 有 效 地 维持 图 像 纹理 细节 ， 文 献 [5] 提 出 观 评价 数据 ,而 且 纹理 清晰 边缘 平滑 , 具有 较 好 的 视觉 效果 。 
有 理 图 像 插值 算法 ， 该 算法 通过 优化 参数 拟 合 待 插值 曲面 ， 

实现 了 纹理 细节 的 保持 及 任意 比例 的 图 像 放 大 。 基 于 文献 [5] 。 1 ” 双 变量 有 理 插值 函数 

提出 的 有 理 插值 函数 ,文献 6-8] 分 别 对 其 进行 了 改进 ,改进 后 给 定 一 幅 大 小 为 mx 的 图 像 ， 即 原始 图 像 区 域 为 
的 插值 算法 能 有 效 地 维持 重建 图 像 的 纹理 细节 ， 但 对 于 图 像 。 wumslal 。 假设 将 图 像 模 向 放大 到 宽度 为 及， 纵向 放大 到 高 
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度 为 y»， 即 放大 后 的 图 像 区 域 为 W:[L.;1,y,]。 令 4,j) 表示 原 
始 像素 点 的 坐标 ， 万 表示 该 点 对 应 像素 值 。(%,3y)) 表示 该 点 
在 放大 后 对 应 的 坐标 , 其 中 dm 人 ns,1< yyj<»)。dijy 和 
ej) 为 点 人) 处 的 偏 导数 Gy)/& 和 (x,3))/60y; 。 令 
对 于 待 插值 区 域 上 的 任意 待 插值 
点 7)s[Dayj yl， 令 0=(-x) 记 7=(07-y)/ 。 对 每 一 
个 y=yj,j=42,.…n， 创 建 * 方 向 的 插值 曲线 ， 表 达 式 为 


已 CD ， 


付 鹏 试 ， 等 : 


l= yn— yo 


hi = Xn 一 故 ， 


(和 DE 1 
Tg 
其 中 : 
FWD=1-0Pf; +01-0° VD+OU- OW +O fi 
qj(X)=(1 -0 +0(-00;+0° 
Vi (2) = (Qi 学 Df; + hd +0(— (1-0a,; finj = —hdi,) 
Wi (xX)= (Qi 号 Dfin,j + hd +0(—(1— Oa; (fin = fj 过 hdinj) 
Qij >0 。 
式 () 定 义 的 插值 函 数 BC5) 满足 BC5) 至 fii » P(xXin) = fi » 
PP’ (%)= dij 9 PP’ (xi) = din,j o 
令 A3=(fj -fi)/h，A=(fim-fij)/1i 。 定 义 
OY ea 008m 
ht+h 
dj = A -AN -AD) O) 
ht+h 
二 A - hn (7) AAW 
diy = 人 ma (An — Am-2.)) 
对 于 每 一 个 点 对 让 i=1,2,.m,j=1,2,.n ,使 用 x 方向 的 函 


数 Pr (x) 来 定义 [Xs Xi; yj, yj] 上 的 函数 Pj(x,y) ， 表达 式 如 下 : 


fj( , 列 
CD) ，_ @) 


i127 7s1,2,..n 
和 07) 


Pj(x,y)= 


其 中 : 
tj(x,y)= -np D+n -7 Vt HW +m P(X) 
gi))= 1-7) +n71 -mp + 
Vj=(B; +DP' +LG (XK) +B (x, y) 
W.;=(B.; + DP (XN) -Lp (A) + Wij (x,y) 
G(X)=(1-00+40+90°)e, +0°(6—80+30)e,,s=j,j+1 
hi = -nd -DB + PAW -P(X -L197)) 
yi P= -nd -npB + PD-—Pin(W) -Lp(W) Bj;>0 。 
式 (3) 定 义 的 函数 B(x,y) 为 双 变 量 有 理 插值 函数 ， 满 足 
Pj ys)= fy) , OP%,y)/x=d,, , OP (x,,))/0y=e,, ， 
r=i,itl;s=]J,j+l1,。 


该 插值 函数 表达 式 由 插值 数据 和 形状 参数 共同 确定 。 
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结构 
部 分 
构造 


所 示 。 首 先 提取 低 分 辩 率 图 
统计 法 设 定 自 适 应 阔 值 将 图 
部 分 ; 然后 基于 边缘 数据 的 
E 离 ， 初 步 确定 待 
插值 点 水 平和 垂直 方向 梯度 角 ， 据 此 计算 已 知 插值 点 
待 插值 点 水 平方 向 梯度 权重 和 垂直 方向 的 梯度 权重 。 
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纹理 细节 和 边缘 结构 保持 的 图 
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作为 图 像 边缘 的 两 个 重要 特性 ， 他 们 在 插值 过 程 中 边缘 
的 保持 上 起 到 了 重要 的 作用 。 因 此 ， 本 文通 过 图 像 边 缘 
的 局 部 不 对 称 性 和 梯度 信息 调整 像素 点 的 空间 距离 使 得 
的 插值 曲面 更 接近 于 理 且 。 

纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算 法 流程 图 如 图 1 


值 点 空间 


的 八方 向 梯度 特征 ， 采 用 简 
x 域 划 分 为 边缘 和 非 边 缘 两 
部 不 对 称 性 调整 边缘 区 域 待 
值 点 坐标 。 在 此 基础 上 分 别 


于 沪 了 芒 条 浇 豆 


区 i 
En 


根 


让 半天 涝 填 寺 二 


获得 的 水 平方 向 梯度 权重 和 垂直 方向 的 梯度 权重 实现 待 插 
点 自 适 应 梯度 的 空间 距离 修正 ， 获 得 边缘 区 域 待 插值 点 修 
后 的 坐标 ， 最 后 基于 边缘 和 非 边缘 部 分 的 结构 特征 ， 在 图 像 
的 非 边缘 部 分 采用 多 项 式 插值 模型 进行 插值 ， 边 缘 部 分 将 修 
下 后 的 待 插值 点 坐标 取代 原 坐 标 代入 有 理 函数 插值 模型 实现 
边缘 区 域 插值 。 
基于 数据 局 部 不 对 
称 性 调整 空间 距离 
梯度 调整 获取 水 平 
和 垂直 方向 梯度 角 
梯度 方向 权重 的 空 
间距 离 调 整 
否 
图 1 纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算 法 流程 
Fig. 1 Gradient adaptive rational function interpolation flowchart 
2.1 图 像 区 域 划分 
图 像 的 边缘 部 分 和 非 边缘 部 分 具有 不 同 的 结构 特征 。 边 


缘 部 分 结构 复杂 , 像素 值 变动 较 大 , 而 非 边缘 部 分 结构 简单 ， 


区 状 参数 的 取 值 不 同时 ， 该 函数 的 表达 形式 也 不 相同 。 
Qj=B;=1， 则 该 函数 为 多 项 式 插值 函数 。 


2 ”纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算法 


其 于 双 变 量 有 理 插值 函数 的 图 像 插值 算法 能 有 效 地 保持 
重建 图 像 的 纹理 细节 信息 的 特点 ， 本 文 借鉴 双 变 量 有 理 插 值 


蚤 数 ， 构 建 双 变 量 有 理 插值 模型 。 但 是 ， 该 双 变量 有 理 函 数 


区 域 葛 
scharr 


过 简单 统 i 


像素 值 变 动 较 小 。 为 提高 图 像 插 值 的 精确 度 ， 不 同 的 部 分 需 
要 采用 不 同 的 插值 模型 。 


本 文通 过 基于 scharr 算 子 的 八方 向 边缘 检测 算法 将 图 像 


| 分 成 边缘 部 分 和 非 边缘 部 分 。 首 先 ， 利 用 八方 向 的 


模板 与 图 像 做 平 


看 卷 积 提 取 图 像 的 梯度 信息 。 随 后 通 


| 法 计算 图 像 每 一 像素 点 对 应 的 平均 像素 值 作为 该 


点 的 冰 值 ， 计 算 公式 如 下 。 最 后 比较 图 像 每 点 对 应 梯度 值 和 


卫 

插值 模型 在 方向 上 是 非 对 称 的 ， 如 果 直 接 将 该 插值 模型 作用 
于 图 像 插值 ， 插 值 结果 的 非 边缘 区 域 插值 效果 较 好 ， 而 边缘 
区 域 会 产生 饮 齿 现象 。 
巧 于 双 变 量 有 理 插 
重建 图 像 的 纹理 细节 信 


值 函 数 的 图 
息 的 特点 ， 


像 插值 算法 能 有 效 地 保持 
本 文 借鉴 双 变量 有 理 插 值 


区 域 


函数 ， 构 建 双 变量 有 理 插值 模型 。 但 是 ， 该 双 变量 有 理 函 数 
插值 模型 在 方向 上 是 非 对 称 的 ， 如 果 直 接 将 该 插值 模型 作用 
于 图 像 插值 ， 插 值 结果 的 非 边 缘 区 域 插值 效果 较 好 ， 而 
区 域 会 产生 饮 具 现象。 

自然 图 像 中 ， 图 像 边 缘 数 据 由 于 获取 过 程 中 的 低 通 滤波 
旺 sigmoid 曲线 分 布 ， 具 有 局 部 不 对 称 的 特性 [0。 多 


其 中 : 


o 


T 


ij 


ye 


e(x,y)=max{|e. ,le, 1} ; 


阔 值 的 大 小 关系 ， 判 断 该 点 所 属 部 分 为 边缘 区 域 还 是 非 边缘 


> > esy)7Coy) 


,G-2<x<i+2.1-2<y<j+2) 


(4) 


e@=f(x-Ly)-fx+thy); 6 =fXy-D- 


Jooy+D ，f(wy) 是 (x,y) 点 的 像素 值 ,e, 为 水 平方 向 的 灰 度 差 
值 ,。, 为 垂直 方向 的 灰 度 差 值 ，e(x,y) 为 水 平和 垂直 方向 中 最 


大 的 灰 度 差 值 。 
应 阀 值 的 八方 向 scharr 边缘 检测 算法 能 够 检测 出 不 


自 适 


同 角度 的 线段 ， 检 测 到 的 边缘 部 分 不 仅 具 有 较 好 的 连续 性 ， 


度 则 体现 了 强度 的 变化 ， 是 进行 未 知 信息 估计 的 重要 依据 2 


而 且 轮 廓 较为 完整 。 图 2 为 lena 采 上 


该 方法 检测 边缘 的 结果 ， 
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可 以 看 出 ， 该 算法 能 够 精确 地 检测 到 图 像 的 有 效 边缘 。 (frpfmim) 分 别 代 入 式 (5) 获 取 待 插值 点 (i, 站 在 一 维 垂直 方 
| 向 上 修正 的 空间 距离 sn 、sj; 。 图 4(c) 显 示 了 (i, 让 在 水 平和 


垂直 方向 上 进行 修正 获得 的 空间 距离 。 


/人 fa+L)D 
31 


JJ]) fli+l,)) f0,N) JrG+L 让 
SG 一 51 一 


pa 2 图 像 边 缘 检 测 结果 JGTJ+D JG+L7T+D 


(a) 水 平方 向 空间 距离 调整 


Fig.2 Image edge detection result 


2.2 自 适 应 梯度 的 空间 距离 调整 


fl f+1,7) 


依据 局 部 不 对 称 性 和 梯度 特征 在 图 像 边缘 结构 保持 中 的 TO St eh eli 
重要 作用 ， 本 节 提 出 将 基于 局 部 不 对 称 性 与 梯度 特征 相 结 合 ; Ce 
的 空间 距离 修正 ， 调 整 待 插值 点 位 置 的 方法 。 FOI+D) FO+hitD 
2.2.1 基于 局 部 不 对 称 性 的 空间 距离 修正 SE 

于 低 通 滤波 边缘 数据 呈 sigmoid 曲线 分 布 ， 因 而 可 用 图 4 二 维 空间 距离 调整 

图 3 所 示 sigmoid 曲线 作为 图 像 的 边缘 模型 。 在 :和 ji+l 点 对 Fig.4 Two-dimensional space distance adjustment 
应 像素 值 分 别 为 fi 和 f,， 而 竺 插值 点 i 的 像素 值 需要 通过 插 最 后 将 一 维 水 平方 向 和 一 维 垂直 方向 上 对 待 插 值 点 进行 
值 算法 计算 得 来 。 一 般 插值 算法 在 计算 竺 插值 点 i 的 像素 值 调整 的 空间 距离 整合 到 二 维 空间 中 。 如 图 5 所 示 ， 根 据 修 正 


时 直接 将 待 插值 点 坐标 代入 插值 函数 中 ， 仅 仅 是 考虑 到 i 与 。” 后 (四 一 维 水 平方 向 上 的 距离 s+ 、ss ， 和 一 维 垂直 方向 上 的 
己 知 点 i 之 间 的 距离 s， 忽略 了 其 他 因素 。 致使 待 插 值 点 的 实 。 距离 s,s 、s; 得 到 如 下 几何 关系 , 整理 后 得 到 式 (7), 即将 si、 


际 位 置 i 与 理想 位 置 i(real) 产生 距离 偏差 ， 得 到 的 像素 值 / 5 、 5 sa 代入 式 (7), 获得 二 维 空间 上 修正 的 待 插值 点 Gi, 让 
与 理想 像素 值 fen 不 符 ， 插 值 结果 不 准确 。 的 水 平 距 离 % 和 垂直 距离 9” 。 
人 图 (5) 得 到 如 下 几何 关系 : 
剖 “ 晤 一 而 
1 So 一 5 
5) _ SS 
下 Sj2 — Sj 
求解 得 : 
图 3 图 像 边缘 数据 sigmoid 曲线 分 布 ， Sa + (Si2 一 Sa)si Sn + (Sj2 — Sn)si 0) 
Fig.3 Sigmoid curve distribution of image edge data 1— (sj — S$)(Si2 — $a) 1—(sj2 — $n)(Si2 — $n) 


为 得 到 边缘 部 分 更 准确 的 插值 结果 ， 需 要 调整 待 插值 点 
与 已 知 点 之 间 的 距离 。 如 图 3 中 ， 当 待 插值 点 i 与 插值 点 i 之 
间距 离 为 时 对 应 像素 值 为 ,而 当 调 整 i 与 i 之 间 的 距离 为 
s 时 对 应 像素 值 为 fj， 即 为 理想 像素 值 。 
基于 图 像 的 局 部 数据 不 对 称 性 09 对 待 插值 点 与 插值 点 
之 间 的 空间 距离 进行 调整 。 以 一 维 方向 上 待 插值 点 的 空间 距 
离 调 整 为 例 ， 设 为 待 插值 点 像素 值 ， 其 坐标 为 xi， 相 邻 点 


的 化 标 分 别 为 i 左 )，i+1( 右 )。 定 义 和 i 之 间 的 距离 为 : 
s=i -i， 将 s 代 入 式 (5) 计 算得 修正 后 的 空间 距离 s 。 图 5 一 维修 正 距 离 转 换 向 二 维 空间 距离 转换 

$s =s—kAs(s—1) (5) Fig. 5 One-dimensional modified distance transformation to 
其 中 ，A 由 式 (6) 定 义 ， 表 示 待 插值 点 邻 域 的 几何 相似 性 。 若 two-dimensional spatial distance transformation 
A=0 则 表示 左右 两 边 像素 对 称 ; 若 A>0 则 表示 竺 插值 点 7" 的。 2.2.2 获取 图 像 梯度 角 
位 置 更 接近 右边 像素 ， 因 而 右边 像素 应 该 具有 较 大 的 加 权 系 设 原始 图 像 为 低 分 辨 率 图 像 ， 首 先 通过 双 三 次 插值 获取 
数 。A<0 时 表示 竺 插值 点 的 位 置 更 接近 左边 像素 ， 左 边 像 ” 初始 的 高 分 辩 率 图 像 ， 其 次 采用 sobel 算 子 分 别 计算 初始 高 


素 应 该 具有 更 大 的 加 权 系 数 。L 在 这 里 表示 灰 度 级 的 最 大 值 ， 分 辩 率 图 像 在 水 平方 向 和 垂直 方向 的 梯度 和 梯度 角 ， 随 后 根 
8 位 图 像 中 取 256， 以 保证 4e[-1,1] 。Kk 为 修正 因子 ， 通 常 取 ”” 据 中 心 像素 梯度 角 对 初始 高 分 辨 率 图 像 梯度 调整 ， 得 到 分 布 
k=1 或 2。 趋 于 合理 的 高 分 辨 率 图 像 在 水 平方 向 的 梯度 G 及 垂直 方向 
太太 -La 下 (0) 的 梯度 G6; 。 最 后 计算 水 平方 向 所 成 的 梯度 角 wD ， 其 中 
L-1 9(i)=arctan|Gi1(G))| ， 及 垂直 方向 所 成 的 梯度 角 90) ， 其 中 
如 图 4 所 示 ， 将 一 维 方向 上 的 空间 距离 调整 扩展 到 二 维 9()=arctan|G) /(G)| 。 
空间 中 。 首 先 在 一 维 方 向 调整 待 插值 点 (i, 站 的 空间 距离 。 斩 以 插值 图 像 中 任意 一 像素 点 Gi)) 在 水 平方 向 的 梯度 调整 
图 4(a) 所 示 ， 获 取 待 插值 点 在 一 维 水 平方 向 上 修正 的 空间 距 ”为 例 ， 说 明 本 文 所 用 梯度 调整 方法 。 在 以 GD) 点 为 中 心 、 大 
离 。 依 据 水 平方 向 上 的 像素 点 (让 ，(i+l, 四 及 其 对 应 像素 值 。 小 为 5x5 的 邻 域 窗口 Nj 内 ， 记 邻 域 窗口 Nj 内 的 任意 一 点 像 
i; ， finj 代入 式 (5) 获 得 待 插 值 点 水 平方 向 上 的 修正 距离 s,。 ” 素 点 (x,y) 水 平方 向 的 梯度 与 (i,)) 在 水 平方 向 梯度 之 间 的 夹 
同 理 依 据 水 平方 向 上 的 像素 点 Gj+D，(i+tbj+D 及 其 对 应 像 ” 角 为 B,， 找 到 BB, 小 于 45 的 所 有 领域 像素 ， 共 计 个 ， 由 公 
素 值 上 ma ，finn 代入 式 (5) 获 得 待 插值 点 水 平方 向 上 男 一 个 修 。 式 (8) 计 算出 这 k 个 像素 的 水 平 梯度 均值 6,， 并 将 (点 的 
正 距离 % 。 随 后 如 图 4(b) 所 示 ， 通 过 自家 方向 上 的 两 组 像素 水 平 梯度 值 用 Gs 来 代 蔡 。 同 理 (ij) 点 的 垂直 方向 梯度 通过 计 
点 (GG 访 (6j+D)) ， (Gt4, 记 ,(i+1,j+D)) 及 其 对 应 像素 值 (fij,fijn) 和  ” 算 后 用 G; 代 蔡 。 
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0 @ 
其 中 Gx, 表示 (%,y) 点 对 应 水 平方 向 上 的 梯度 值 。 
2.2.3 梯度 自 适应 的 空间 距离 调整 
根据 数据 局 部 不 对 称 性 进行 空间 距离 调整 ， 获 得 新 的 待 
插值 点 坐标 5 和 s;， 只 考虑 了 边缘 数据 信号 分 布 对 空间 距离 
的 影响 ， 未 能 充分 地 挖掘 图 像 中 所 含 信息 。 由 于 图 像 边缘 处 
像素 的 梯度 特征 非常 显著 ， 因 此 图 像 的 梯度 信息 对 于 重建 图 
像 边缘 结构 的 保持 有 着 重要 作用 。 本 文 算法 在 基于 局 部 不 对 
称 性 的 空间 距离 修正 后 结合 梯度 特征 ， 再 次 调整 图 像 的 空间 
距离 ， 获 得 新 的 插值 点 坐标 。 
观察 分 析 得 ， 邻 域 像素 点 对 应 像素 值 沿 中 心 点 像素 的 
梯度 方向 变化 最 为 明显 。 即 在 中 心 像 素 点 的 邻 域内 ， 邻 域 像 
素 点 像素 值 在 中 心 像素 点 的 梯度 方向 上 的 投影 距离 远近 ， 表 
明了 两 者 的 空间 相关 性 强 弱 。 基 于 此 ， 可 通过 邻 域 像 素 点 在 
待 插值 点 梯度 方向 上 投影 距离 的 远近 ， 确 定 待 插值 点 坐标 的 
偏 移 程 度 ， 调 整 待 插值 点 与 邻 域 像素 点 之 间 的 距离 ， 以 获得 
符合 插值 曲面 分 布 的 待 插值 点 坐标 。 本 文 算法 首先 分 别 计算 
邻 域 像素 点 在 待 插值 点 水 平 梯度 方向 和 垂直 梯度 方向 上 的 投 
影 距 离 ， 其 次 依据 该 投影 距离 获取 邻 域 像素 点 在 待 插值 点 水 
平 梯度 和 垂直 梯度 方向 上 的 梯度 特征 权重 ， 再 次 依据 水 平和 
垂直 方向 的 梯度 特征 权重 计算 待 插值 点 在 水 平和 垂直 方向 上 
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代理 


| 


纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算法 
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度 方向 上 的 


域 像 素 点 在 待 插值 点 水 平 梯度 方向 上 的 投影 距离 ds， ds ， 
da ，du， 代入 权 值 函数 公式 (11) 分 别 计算 邻 域 像素 点 在 待 插 


值 点 水 平方 


其 中 “为 党 


于 归 一 化 


Fig. 6 


再 次 ， 


句 邻 域 像素 


向 邻 域 像素 点 的 偏 移 距 离 ， 最 后 待 插 值 点 在 水 平方 向 和 垂直 
方向 上 的 偏 移 距离 拟 合 到 二 维 空间 ， 调 整 待 插值 点 位 置 ， 获 
得 自 适 应 梯度 修正 后 的 待 插值 点 坐标 (2*,j”*) 用 以 计算 待 
插值 点 (i,7 对 应 像素 值 。 
以 水 平方 向 上 待 插值 点 空间 距离 修正 为 例 进 行 说 明 。 首 
先 , 计算 邻 域 像素 点 在 待 插值 点 水 平 梯度 方向 上 的 投影 距离 。 


度 系 数 Un ， 
方向 上 与 令 


最 后 ， 
空间 ， 调 整 


如 图 6 所 示 ， 根 据 待 插值 点 垂直 方向 华 标 六， 构造 一 个 1 
六 Jo 站 的 所 有 点 构成 的 平面 1。 将 邻 域 像素 点 对 应 像素 
值 垂直 投影 到 该 平面 上 。 由 待 插值 点 水 平 华 标 i， 双 三 次 插 


梯度 方向 的 
值 点 空间 距 


值 获取 的 初始 像素 值 f0,j) 及 2.2.2 节 获 取水 平方 向 梯度 角 
9(i) ， 在 平面 1 上 构造 一 条 过 点 (i,f(i, 门 ) 且 斜率 为 tan(8()) 的 
直线 。 如 下 所 示 为 该 直线 的 函数 表达 式 


的 待 插 值 点 


4 表示 邻 域 像素 点 投影 到 待 插值 点 梯度 方向 上 的 距离 。m 用 


空间 修正 距 


梯度 特征 权重 。 定 义 权 值 函数 如 下 所 示 。 根 据 邻 


向 上 的 梯度 所 占 权 重 Un, U2, Ur, Uo 


u= (11) 


m 


数 , 用 于 控制 指数 函数 的 衰减 ,一 般 设置 值 为 0.2。 


dx3 drivy{ (itl, yj) 
dx 和 4 


f(i, j+1) 


f(i+1 计 1) 


图 6 待 插值 点 水 平 梯度 方向 投影 距离 


Projection distance in horizontal gradient direction of 


points to be interpolated 
所 得 梯度 特征 权重 计算 待 插值 点 在 水 平方 向 上 
点 的 偏 移 距离 。 利 用 待 插值 点 水 平方 向 上 局 部 梯 
ws ，ui，oa 代入 式 (12)， 获 得 待 插值 点 一 维 水 平 
或 像素 点 @7 的 距离 H， 及 (+1 四 的 距离 $5。 


站 Un(s— 58) +uns: ” (12) 
将 待 插值 点 在 一 维 方 向 上 修正 的 距离 拟 合 到 二 维 
待 插值 点 的 水 平 坐标 。 根 据 获得 的 待 插值 点 水 平 
维 空 间 修正 距离 s+ ，s; 和 通过 垂直 方向 上 待 插 
离 修正 计算 获得 的 待 插值 点 垂直 梯度 方向 的 一 维 
离 sn ， 5s 通过 式 (7) 将 其 转换 为 二 维 空间 上 修正 
水 平 距 离 ss。 这样 完 成 了 待 插值 点 在 水 平方 向 上 


US; W125 


$i2 二 
U2 (5 — $i) + U5; 


自 适应 梯度 


坐标 real 。 


的 空间 距离 修正 ， 获 得 了 修正 后 待 插值 点 的 水 平 


Jo 门 =E+D 2.3 ”基于 有 理 函 数 的 图 像 插值 
k=tan(g(i)) (9) 节 构 建 双 变 量 有 理 插值 模型 实现 图 像 插值 ， 该 模型 通 
Ca CE 过 参数 调节 可 在 有 理 函 数 模型 和 多 项 式 模 型 之 间 转 换 。 有 理 
其 中 x 为 该 直线 上 任意 一 点 的 水 平 坐标 ， 了 (x, 站 为 该 直线 上 函数 模型 能 够 在 很 小 的 区 域 使 得 图 像 曲 面 达到 极 值 ， 有 利于 
水 平 坐 为 x 的 点 对 应 的 像素 值 。 表 示 该 直线 的 斜率 , pC) 为 ”处 理 像 素 值 波动 较 大 的 区 域 。 而 多 项 式 模型 能 够 将 插值 曲面 
待 插 值 点 ,7 对 应 的 梯度 角 。4。 为 该 直线 的 截 距 。 控制 在 一 定 范围 内 变动 ,有 利于 处 理 数 据 变动 比较 小 的 区 域 。 
直线 的 函数 表达 式 ， 计 算 该 平面 上 邻 域 像素 点 ”考虑 图 像 边缘 部 分 和 非 边缘 部 分 的 结构 特征 ， 图 像 边缘 区 域 
GD) tl f+) ，(i fj+D)) ，(i+1,f(i+1,j+)) 在 待 插 ” 采用 有 理 模型 ， 而 非 边 缘 区 域 采 用 多 项 式 模型 。 
值 点 水 平 梯度 方向 上 的 投影 距离 ds ，d., ，ds ，dws 。 如 图 6(b) 如 图 7 所 示 ， 假 设 4 个 黑色 像素 点 所 围 成 区 域 1 为 待 插 
所 示 ， 为 该 水 平 切面 ( 即 平面 1) 的 平面 示意 图 。 以 平面 1 上 域 ” 值 区 域 ,4 个 红色 像素 点 为 待 插值 点 。 由 式 (3) 可 知 ， 根 据 待 
像素 点 ,fi 四 在 竺 插值 点 水 平 梯度 方向 上 的 投影 距离 d 的 ” 插值 区 域 1 周围 的 12 个 像素 点 可 以 构造 出 7 的 插值 曲面 。 位 
计算 过 程 为 例 进行 说 明 。 首 先 依据 直线 的 函数 表达 式 (9), 计 ”于 了 中 的 任意 一 点 均 可 由 7 所 在 插值 曲面 的 函数 表达 式 表示 。 
算 邻 域 像 素 点 Gf67) 到 直线 的 距 离 。 若 1 所 在 区 域 为 边缘 部 分 将 经 空间 距离 调整 后 的 待 插值 点 
dis_1 直 kxi- (站 +b|/NR+1 ; 随后 计算 该 点 与 待 持 值 点 之 间 的 。 ”坐标 代入 所 属 插 值 曲面 函数 表达 式 , 获得 该 点 对 应 的 像素 值 。 
距离 dis_2=VG-D+0(, 站 -J(j) ; 最 后 获取 邻 域 像素 点 。 反之， 为 非 边缘 区 域 时 ， 令 式 G3) 中 wy 和; 均 为 1， 上 述 模 
(fi 四 到 待 插值 点 的 投影 距离 4 =Vas -2 -ds_P 。 整 理 得 平 。 型 转换 为 多 项 式 模 型 。 直 接 代入 待 插值 点 坐标 ， 获 得 该 点 对 


而 上 邻 域 像 素 点 到 待 插值 点 的 投影 距离 计算 公式 为 式 (10)。 
分 别 将 平面 上 其 他 邻 域 像素 点 G+Lf 了 Ot) ，(i,f(i,j+D))， 
G+lfG+7+D) 代入 式 (10) 计 算 其 在 待 插值 点 水 平 梯度 方向 
上 的 投影 距离 Co ，ds ， dw 。 


4=WJG-PD+OG 门 -GExi- GD+bANG+D2 (10) 
其 中 为 邻 域 像 素 点 的 水 平 坐标 ， ,7 为 邻 域 像素 点 对 应 像 


素 值 。 
其 次 ， 依 据 投影 距离 获取 邻 域 像素 点 在 待 插值 点 水 平 梯 


应 的 像素 值 
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Fig.7 Image interpolation 
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吉 果 分 析 


8 所 示 , 本 文选 择 8 
隔行 隔 列 的 方法 对 测试 图 
重建 的 实验 结果 ， 选 择 


pa 


图 像 作为 测试 
像 进 行 下 采样 
有 代表 性 的 bicubicD， 


， 针 对 


本 文 医 
NEDII， 


像 进行 实验 。 


区 
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Lee” sl!5] 


， 梯 度 优 化 的 有 到 
值 算 法 Gational function interpolation algorithm 
F NSCT 的 
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区 域 


适应 图 像 插值 算法 (an region adaptive image interpolation 


algorithm based on the NSCT,NSCTRAD09 分 别 从 客观 数据 、 


主观 效果 时 间 复 杂 度 三 个 方 


和 i 对 重建 效果 进行 比较 。 


表 1 


构 相 


以 性 (structural similarity index， 
PSNR 是 通过 衡量 标准 


果 2 


间 的 像素 差异 , 对 医 


标准 
像 


对 


区 


视觉 感知 的 质 


二 


量 评价 
NEDI, DFDF, RASI, NARM, Lee’ 


SSIMI7J 是 


图 像 平均 质量 的 评价 ， 


像 


根据 表 1 的 实验 


展示 了 不 同 插值 算法 峰值 信 品 比 (peak signal to noise 
ratio，PSNR) 和 结 
的 实验 对 比 结果 。 其 中 


SSIM) 


图 像 和 插值 结 


根据 


像 与 实验 结果 间 结 构 、 亮 度 和 对 比 度 的 相似 程度 ， 对 
结果 得 与 bicubic 
s，FDI，GORI 方法 相 
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有 最 高 的 PSNR 均值 和 SSIM 均值 
而 在 客观 数据 上 本 文 算法 优 人 


表 1 不 同 插值 算法 的 PSNR 和 SSIM 对 比 
Tab.1 Comparison of different interpolation algorithms on PSNR and SSIM 
% Baboon Fence Barbara Girl Wall Papper Sky Airplane Average 
i psnr SSIM psnr SSIM psnr SSIM psnr SSIM psnr SSIM psnr SSIM psnr SSIM psnr SSIM psnr SSIM 
Bicubic 22.66 0.84 21.82 0.73 23.39 0.90 29.40 0.95 24.27 0.87 29.77 0.95 28.42 0.85 25.78 0.90 25.56 0.87 
NEDI 22.56 0.87 22.79 0.79 22.36 0.85 30.19 0.98 23.94 0.88 30.67 0.97 28.87 0.92 26.95 0.91 25.65 0.89 
DFDF 22.81 0.86 23.25 0.80 24.64 0.88 30.00 0.96 24.91 0.89 29.69 0.88 30.01 0.91 28.20 0.92 26.67 0.90 
RSAI 23.18 0.87 22.91 0.80 23.37 0.86 30.39 0.97 25.17 0.89 31.91 0.97 29.60 0.91 28.48 0.92 26.62 0.90 
NARM 22.63 0.86 23.08 0.80 23.33 0.87 30.55 0.97 24.63 0.90 30.84 0.96 29.20 0.91 28.29 0.92 26.47 0.90 
Lee’s 22.33 0.86 22.66 0.79 23.03 0.85 30.49 0.97 24.64 0.88 30.85 0.95 29.50 0.91 27.69 0.92 26.25 0.89 
GORI 23.02 0.87 24.83 0.80 24.54 0.88 30.67 0.96 24.87 0.88 31.98 0.96 29.82 0.91 28.31 0.91 27.26 0.90 
NSCTRAI 24.23 0.87 24.92 0.80 23.33 0.88 30.67 0.97 24.90 0.89 33.01 0.96 30.10 0.91 28.48 0.92 27.45 0.90 
本 文 方法 24.30 0.88 25.11 0.81 24.28 0.89 30.89 0.96 24.96 0.89 33.10 0.97 30.71 0.92 28.S2 0.93 27.74 0.90 
图 9、10 表示 不 同 算法 在 针对 图 像 纹理 细节 的 重建 效果 图 9~11 展示 了 不 同 算法 在 图 像 的 纹理 和 边缘 部 分 的 重 
对 比 。 图 9 列 出 了 Barbara 围巾 部 分 经 由 不 同 插值 算法 重建 ” 建 结果 对 比 。 由 上 述 分 析 知 本 文 算法 不 仅 维 持 了 图 像 的 纹理 
的 结果 对 比 。 该 0 有 多 方向 的 纹理 特征 ,NEDI、RAST. 细节 ， 还 保护 了 图 像 的 边缘 结构 。 从 主观 效果 上 看 本 文 算法 
NARM、Lee”s 等 算法 在 该 部 分 图 像 重 建 中 出 现 纹理 断裂 现 ee 建 结 果 
象 ，bicubic、DFDF、NSCTRAI 等 算法 实验 结果 纹理 模糊 ， 机 
并 且 bicubic 实验 结果 锯齿 明显 ，GORI 算法 在 图 像 纹理 区 域 DE i 
重建 效果 较 好 ， 但 是 部 分 边缘 区 域 模糊 。 而 本 文 算法 在 该 图 ee 
像 纹理 部 分 重建 ， 未 出 现 纹理 断裂 和 模糊 等 现象 。 图 10 是 不 Wn 
同 算法 针对 Fence 局 部 重建 效果 的 对 比 ， 由 图 可 知 bicubic、 汐 和 测 
DEFDR NARM 重建 结果 出 现 明显 的 边界 模糊 NEDI RSAT ge Wer ee gz 
Lee”s、NSCTRAI 算法 出 现 不 连续 的 条 纹 ， 同 时 有 边界 失真 图 10 “不同 算法 重建 Fence 对 比 
现象 出 现 ，GORI 算法 重建 结果 在 纹理 清晰 区 域 重 建 效 果 较 Fig. 10 Comparison of Different algorithms reconstruct the Fence 
好 ， 但 对 于 远 处 区 域 重建 结果 模糊 同时 出 现 边 界 模糊 ， 轮 万 
重建 不 清晰 问题 。 而 本 文 算法 不 仅 纹理 特征 被 有 效 的 保留 ， 
同时 ed 没有 出 a 


边界 部 分 重建 ? 


图 9 不 同 算法 重建 barbara 对 比 
Fig.9 Comparison of different algorithms reconstruct the Barbara 
图 11 展示 了 不 同 算法 针对 图 像 边缘 部 分 重建 的 结果 对 
比 。 针 对 Barbara 脸 部 轮廓 边缘 部 分 ，bicubic、NARM、FDI 
算法 重建 图 像 边 缘 具 有 明显 的 锯齿 ，GORI 算法 重建 边缘 纹 
理 扭 曲 ， 在 这 其 中 NEDI 和 本 文 算法 具有 较 平 滑 边 缘 。 
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Fig. 11 
在 实 
es 
间 对 比 。 


际 应 用 


表 2 可 以 看 


但 是 图 像 重 建 效 果 不 理想 ， 轿 


糊 和 失真 ，NEDI 和 Lee” 


不 同 算法 重建 


Girl 对 比 


出 ，bicubic 算法 运行 
算法 和 GORI 次 之 ，NARM 运行 时 间 最 长 。 
综 上 所 述 双 三 次 插值 算法 结构 简 闯 


外 且 


过 程 中 ， 时 间 复 杂 度 是 算法 考量 的 重要 标准 


Comparison of different algorithms reconstruct the Girl 


2 展示 了 不 同 插值 算法 针对 上 述 插值 结 


二 果 的 运 


时 


时 间 复 杂 度 


运行 时 
司 最 短 ， 本 文 


最 低 


理 区 域 容易 产生 模 
s 在 保持 图 像 边缘 结构 方面 效果 突 


录用 定稿 付 有 网坛 ， 等 : 纹理 细节 和 边 
出 ， 但 是 重建 的 纹理 区 域 会 出 现 扭曲 和 变形 ; RASI 算法 在 纹 
理 保 持 方面 同 NEDI 一 样 会 出 现 纹理 扭曲 现象 ， 但 是 区 别 于 


眶 


NEDI， 该 算法 在 图 像 的 边缘 部 分 会 产生 锯齿 ， DEFDF 算法 重 
建 图 像 模 糊 ， 视 觉 效 果 不 理想 ，FDI 重建 结果 会 产生 边界 模 
糊 ，GORI 算法 重 建 结果 纹理 保持 效果 较 好 ， 但 在 图 像 边缘 
部 分 容易 出 边缘 模糊 和 锯齿 等 现象 。 相 比 于 上 述 算法 ， 本 
文 算法 重 妈 边缘 平滑 且 纹 理 清晰 具有 较 好 的 主观 效果 
PSNR 平均 提升 了 0.48dB-2.17dB，SSIM 平均 提升 了 0.004- 
0.028， 获 得 了 较 高 的 客观 评价 数据 ， 且 运行 时 间 较 短 。 
表 2 不 同 插值 算法 运行 时 间 对 比 
Tab.2 Comparison of different interpolation algorithms on run time /s 


算法 


i 


Baboon Fence Barbara Girl Wall Papper Sky Airplane Average 


Bicubic 2 1 3 3 3 2 1 1 2 
NEDI 27 6 25 36 38 24 37 6 25 
DFDF 22 6 25 35 36 20 34 6 23 
RSAI 213 50 215 322 327 198 322 49 214 
NARM 1630 265 1450 22272658 1521 2376 265 1544 
Lee’s 105 22 98 142 151 70 146 20 97 
GORI 5 4 4 5 6 4 5 3 4 

NSCTRAI 30 7 29 40 42 33 40 7 29 
本 文 方法 7 3 7 7 8 4 3 4 6 


4 ”结束 语 


本 文 在 借鉴 双 变量 有 理 函 数 插值 模型 的 基础 上 ， 提 出 了 
一 种 纹理 细节 和 边缘 结构 保持 的 图 像 插 值 算法 。 该 算法 首先 
提取 低 分 辨 率 图 像 的 八方 向 梯度 特征 ， 采 用 简单 A 
自 适应 闹 值 将 图 像 区 域 划 分 为 边缘 和 非 边 缘 两 大 部 分 ， 其 
依据 图 像 边缘 数据 的 局 部 不 对 称 性 和 梯度 特征 调整 边缘 区 域 
待 插值 点 的 空间 距离 ， 获 得 经 空间 距离 修正 后 的 待 插值 点 坐 


水。 最 后 基于 边缘 和 非 边 缘 部 分 的 结构 特征 ， 在 图 像 的 非 边 
缘 部 分 采用 多 项 式 插值 模型 进行 插值 ， 边 缘 部 分 将 修正 后 
待 插值 点 坐标 取代 原 坐 标 代 入 有 理 函 数 插值 模型 实现 边缘 
或 插值 。 本 文 所 提 算 法 不 仅 具 有 有 理 函 数 纹理 细节 保 真 度 
的 优点 ， 同 时 还 充分 利用 了 边缘 数据 局 部 不 对 称 性 和 梯度 
保证 了 图 像 边 缘 结 构 。 实 验 结果 表明 ， 该 算法 不 仅 具 有 


Xl 到 


lt Ad 


较 好 的 视觉 效果 ， 同 时 具有 较 高 的 客观 评价 数据 。 
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